The interpretation of the field angle dependence of the critical current in 

defect-engineered superconductors 
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We apply the vortex path model of critical currents to a comprehensive analysis of contempo- 
rary data on defect-engineered superconductors, showing that it provides a consistent and detailed 
interpretation of the experimental data for a diverse range of materials. We address the question 
of whether electron mass anisotropy plays a role of any consequence in determining the form of 
this data and conclude that it does not. By abandoning this false interpretation of the data, we 
are able to make significant progress in understanding the real origin of the observed behavior. In 
particular, we are able to explain a number of common features in the data including shoulders at 
intermediate angles, a uniform response over a wide angular range and the greater discrimination 
between individual defect populations at higher fields. We also correct several misconceptions in- 
cluding the idea that a peak in the angular dependence of the critical current is a necessary signature 
of strong correlated pinning, and conversely that the existence of such a peak implies the existence 
of correlated pinning aligned to the particular direction. The consistency of the vortex path model 
with the principle of maximum entropy is introduced. 

PACS numbers: 74.25.Sv, 74.25.Wx 



I. INTRODUCTION 

In the field of high temperature superconductor (HTS) 
research and development, a significant hurdle has proved 
to be the anisotropic properties of the available mate- 
rials. For a typical YBa2Cu307 (YBCO) sample, the 
variation in the critical current density with the angle 
of an applied magnetic field, J c (0), can exceed a factor 
ten at high fields [l[ . Successful methods of tackling this 
issue have been developed, most notable of which has 
been the incorporation of artificial pinning centers in the 
form of self-assembled BaZrO,3 nanocolumns Q. These 
have introduced a so-called c-axis peak to the pinning 
profile, enhancing J c for fields applied in the c-direction, 
perpendicular to the film surface. Rather than solving 
the problem, however, this has merely shifted the re- 
gion of concern to intermediate angles, between the two 
peaks, where J c remains low. Adding BaZr03 nanopar- 
ticles creates a more isotropic profile but the effect varies 
substantially with field and temperature Q . In order to 
fully overcome this issue by engineering appropriate de- 
fect structures, we need first to understand the origin of 
the J c (0) dependence. Remarkably, no consensus model 
of J c (0) exists, highlighting the general weakness in es- 
tablished models of J c to reproduce the phenomenology 
of critical currents. 

One proposed origin of at least some of the features 
of J c (0) gaining momentum has been the electronic mass 
anisotropy of the materials [4-6]. However, the physical 
link between this intrinsic material anisotropy and (ex- 
trinsic) properties such as J c has not been made clear. 
While application of this approach to understand the 
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physics of pristine samples of clean superconductors Q 
may have some merit, its direct transfer to explaining 
the pinning in technological materials incorporating ar- 
tificial pinning centers, resulting in "effective electronic 
mass anisotropy ratios" that bear no relation to the fun- 
damental material parameters Q is questionable. We 
therefore begin by reviewing the physics of the electron 
mass anisotropy scaling approach and discuss whether 
this approach is applicable to high- J c pinning-enginccrcd 
samples. 

Additional to this, many recent pinning-optimized 
samples @ exhibit features such as pronounced shoulders 
fioj and fiat angular regions in J c (8) that cannot be 
explained by the naive summation of an afo-plane intrinsic 
pinning peak, a c-axis correlated pinning peak and a ran- 
dom pinning contribution characterized by a geometrical 
anisotropy of whatever origin. The presence of similar 
features throu gho ut the literature from the earliest sam- 
ples produced jl2l [l3j to present-day 'champion' samples 
fl4| means that these can not be dismissed as due to poor 
sample quality or inaccurate measurement, but must be 
acknowledged as genuine effects of an optimized pinning 
landscape. Combine this with the complex temperature 
and magnetic field dependence of J c now starting to be 
revealed flfjj and it is apparent that a deeper understand- 
ing of the effects of a complex flux pinning landscape is 
required. 

Here we demonstrate that an analysis of the statisti- 
cal effects of combining c-axis correlated pinning and ab- 
plane correlated pinning (whether intrinsic, extrinsic, or 
surface in origin) can explain all the features observed 
to date in pinning-optimized samples. Isotropic pin- 
ning resulting from point defects or incorporated second- 
phase nanoparticles contributes to this picture in the 
form of 'bridges' between these defining pinning popu- 
lations. The aim of this paper is to show that the model 



2 



that results from this analysis, termed the vortex path 
model, provides a consistent and detailed explanation of 
the experimental features that have risen to prominence 
in recent studies. Using this model, J c (9) is decomposed 
into statistically distinguishable components that can be 
related back to distinct aspects of the sample microstruc- 
ture that determines J c — a link that is missing in other 
approaches. 

II. MODELS OF J c (6) 
A. Electron mass anisotropy scaling 

The electron mass anisotropy scaling approach has re- 
cently become a common method of analyzing J c {9) data 
obtained from diverse species of superconducting mate- 
rials [l6l - fl9| . In some cases misinterpretation of the pro- 
cess has lead to simple fits of the Ginzburg-Landau mass 
anisotropy function to J c (9) being employed [20| . It is 
therefore worth reviewing the reasoning behind this ap- 
proach to understand the limits of its applicability. If the 
scaling approach is valid then this affects how we should 
interpret the various features of J, (8). 

The sca ling approach was originally proposed by Blat- 
ter et al. [211(a)] . with a similar formulation being inde- 
pendently proposed by Hao and Clem [22| . The idea is to 
apply results derived for an isotropic superconductor to 
an anisotropic superconductor through a rescaling of the 
magnetic field, temperature, characteristic lengths, etc. 
The scaling rules are obtained by finding a transforma- 
tion which recasts the Gibbs free energy of an anisotropic 
superconductor into the form for an isotropic material. 
The primary result is then a scaling rule such as 

Q(9, H, T, e, X L ,e, S) = s Q Q(e e H, T/e, £, X L , 8/e). (1) 

Blatter et al. then apply this rule to calculate criti- 
cal currents for the case of single-vortex weak collec- 
tive pinning, that is pinning by randomly distributed 
point pins. The results arc jj = J^, independent of 
8, for the in-plane critical current density (as usually 
measured in thin film experiments), where J£ is the in- 
plane J c with the field applied in the c-dircction, and 
= egJ£ for the out-of-plane critical current density, 
with Eg = e(9) = cos 2 9 + e 2 sin 2 8, where e 2 = m a b/m c is 
the electronic mass anisotropy. For consistency through- 
out our paper, and in contrast to Blatter, we have taken 
6 = to be perpendicular to the afe-planes. 

The result jj = constant for low fields is consistently 
ignored by those who use the scaling approach to anal- 
yse J c (9) data. We are not aware of any report claim- 
ing to have confirmed this result experimentally, thereby 
validating the scaling hypothesis. Many samples re- 
ported in this paper and elsewhere have wide regions of 
jj = constant at low field. However, these samples also 
have correlated c-axis defects such as grain boundaries, 
threading dislocations and twin plane intersections and it 



would be puzzling to assume these made no contribution 
to pinning in these circumstances. 

The scaling procedure that is widely used in the litera- 
ture is based on an expression postulated in a subsequent 
review by Blatter et al. [2llfb)] for the high field case, de- 
fined by H{9) > p sb {J c JJ d )H c2 {9) where p sb « 5 and J d 
is the depairing current density. In this case, no physical 
model of J c to which the scaling rules can directly be re- 
lated is proposed. Rather, the authors simply assert that 
at high field, J c will scale in accordance with H — > sgH, 
i.e. 

4(6, H) = J c c (e e H). (2) 

To evaluate the status of this assertion we recall that the 
scaling rule, Eq. ([1]), is a non- unique mathematical trans- 
formation that allows us to write the equilibrium free en- 
ergy for an anisotropic superconductor in the same form 
as for an isotropic superconductor. It is not a physical 
theory where a measurement of a physical quantity may 
actually give a different answer depending on the value of 
another physical variable. In particular, a measurement 
of the field in the sample does not scale with these rules. 
Any change in the critical current must therefore have 
an indirect mechanism, which this approach does not ad- 
dress. The critical current density is a non-equilibrium 
property of the sample; formally, it is a parameter in a 
constitutive equation such as E = Eq(J/ J c ) n . There is 
no formal derivation possible from Eq. ([1]) to Eq. ([2]); 
rather, Eq. @ is a new hypothesis. The justification 
for the adoption of this hypothesis, noted in both [2ll(b])] 
and [22| . is a perceived similarity between the high field 
behavior of J c (9) and the field-scaling rule for the free 
energy. The justification given is empirical, not theoret- 
ical. 

This justification, however, is not well founded. If a 
broad peak in J c (9) is determined by a variation in pin- 
ning density with angle, for example arising from the 
interaction of correlated defects, rather than from the 
mass anisotropy of the material, then we would still ex- 
pect J c (9) data to be approximately scaleable in this 
way. Typically, for HTS samples, critical currents in 
large out-of-planc fields behave as J c (H\\c) cx H~ a . If, 
for an in-plane field, we have a higher density of cor- 
related pinning defects but otherwise similar physics we 
would expect J c (H\\ab) oc (H/k)~ a , where k is a mea- 
sure of the change in pinning density between the field 
directions. We can scale H — > He(k, 9) so the data co- 
incide at these two angles using an equation of the same 
form as the mass anisotropy, £ 2 (fc, 9) = cos 2 9 + k 2 sin 2 9. 
Further, if J c {9) varies smoothly and broadly then e(k, 9) 
will also give a reasonable fit to all the intermediate angle 
data. The electron mass anisotropy expression can there- 
fore be a reasonable fitting function whether the origin 
of the behavior is in fact electronic mass anisotropy or 
simply pinning density anisotropy arising from correlated 
defects. Experimentalists who employ this approach as- 
sume that correlated pinning can only be effective and 
therefore influence their data over narrow angular ranges. 
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The evidence from samples implanted at different angles 
by heavy ions l23l . |24| or the addition of self-assembled 
nanocolumns |25| suggests this confidence is misplaced. 
Similar to the low field result, those who apply this scal- 
ing also consistently ignore the lower bound to the field 
range over which it is applicable and do not attempt to 
verify any crossover from the low field to the high field 
regime. The lower bound quantified above can be signif- 
icant for HTS since H C 2 is large; at low temperature it 
may easily be tens of teslas, beyond most common ex- 
perimental capability. 

Aside from these criticisms the scaling approach suf- 
fers from one further fatal flaw when applied to pinning- 
optimized samples: it ignores the pinning summation 
problem. If we add random isotropic defects to a sample 
that already possesses a high density of a6-planc pinning 
but much lower c-axis pinning then we expect the effect 
on J c {9) about the c-axis direction to be much stronger 
than the effect about the ab- plane. This is exactly what 
is seen experimentally [H, Hol. l26l| . As individual vortices 
interact with multiple defects, correlated and uncorre- 
cted, the effect on J c of additions to the pinning land- 
scape depends on the pinning already present and the 
scaling approach has no way of taking this into account. 

The use of the high field scaling approach to analyze 
J c {9) has led to a situation where a correspondence be- 
tween the experimental data and the known defect struc- 
tures of HTS films is missing. We think this approach is 
wrong and should be abandoned. In the next section we 
recapitulate a model that can distinguish how different 
defect populations within the sample contribute to J c , 
and that can enumerate how these distinct contributions 
combine to yield the overall J c (9). 



B. The vortex path model 

A statistical mechanical approach to modeling the re- 
sponse of a superconductor to a particular population of 
pinning defects was first proposed in [27} and formally 
presented in 1281. It was subsequently termed the vortex 
path model in [29( since at its core, it considers the multi- 
plicity of possible pinned vortex paths through a sample, 
as determined by its defect structure, to quantify the rela- 
tive pinning strength for fields applied at different angles. 
It has recently been used to describe experimental data 
where the electronic mass anisotropy approach has been 
found inadequate (30j . 

Here we do not repeat the full exposition of the model 
but instead give a general theoretical context as an in- 
troduction. The model is constructed in the spirit of 
a maximum entropy approach [3lj . If we consider any 
macroscopic property of a system, where the macroscopic 
state can be related to a multiplicity of microstates, then 
to predict the behavior of the property we maximize the 
Shannon entropy of the distribution of microstates of the 
system subject to the known constraints. For elemen- 
tary constraints this leads to the ubiquitous functions of 



physics; for example, the Gaussian function is obtained 
upon maximizing the Shannon entropy subject to the 
constraints that the mean and the variance of the dis- 
tribution of microstates are fixed, while the (truncated) 
Lorentzian function is obtained when only the variance is 
fixed [32| . The critical current density of a superconduc- 
tor is just such a macroscopic property arising from the 
multiplicity of microstates in which vortices are pinned 
at a particular angle, and following [28j we can write 



J c (6) = J g(6), 



(3) 



where g(9) is the density distribution of these pinning 
microstates normalized such that J* g(9)d9 = 1 and the 
proportionality factor Jo = J* J c (9)d9. In [28| . a Lorentz 
force definition of the critical current was used to define 
Jo in terms of the pinning force. However, this definition 
of J c is unnecessary, and instead Eq. [3] itself serves to 
define the critical current. We can then maximize the 
Shannon entropy 



H 



g(6)Hg(6))d6 



(4) 



with appropriate constraints in order to predict J c {9). 
Without any constraints the result is a uniform distribu- 
tion and we obtain a constant J c , 



MO) = i 

7T 



(5a) 



For a sample containing correlated defects within the 
afe-plane and along the c-axis we know that J c (9) may 
exhibit peaks centered on the directions of the defects. 
These peaks are broadened by the vortices' interactions 
with the orthogonal defect populations. The vortex path 
model introduces constraints on g(9) by requiring the 
vortex microstates to follow a path which interacts with 
these defects. In the usual geometry for J c measurements 
we define 9 = tan (y/z), where y lies within the plane 
of the sample and z is perpendicular to it. For a peak 
centered on y, we constrain the pinned vortex path to 
comprise m steps of average y-distance A, so y = mX is 
the vortex length in the y-direction, and m steps of length 
z% in the z-dircction, where the distribution of lengths 
z% is p(zi). The total z displacement is then given by 
z = ^2zi, with a distribution p(z). The principle of max- 
imum entropy can be applied and we hypothesize entropy 
maximizing distributions ofp(z) and hence g{9), leading 
to different results for J c (9) through a transformation of 
random variables from (y, z) to 9. If the distribution p(z) 
is Gaussian, we obtain 



M9) 



Jo 



27rr sin 2 9 



exp 



f 



2r 2 tan 2 9 J ' 



while if p{z) is Lorentzian, we obtain 



f 



Jc(9) = - 



J r 



7r cos 2 9 + T 2 sin 2 9 



(5b) 



(5c) 
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We term Eqs. (|5b[l and ((5c)) angular Gaussian and an- 
gular Lorentzian functions, respectively. In Eq. (|5b[) . 
r = a/y/mX, where ma 2 is the variance of the Gaus- 
sian, while in Eq. (|5c[) , V = 7/ A, where 7717 is the scale 
factor of the Lorentzian. Note that Eq. (|5c|) incorporates 
the special case T = 1 for which it reduces to the uniform 
function of Eq. (f5"aj) . 

To explain this further, the choice of a Gaussian dis- 
tribution for p(z) indicates that there exists a popula- 
tion of orthogonal defects that broaden the peak with 
the only added constraints that p{z) has a mean and a 
variance, i.e. that the lengths of pinned vortex segments 
in the z-direction form a normal distribution. If multiple 
defect populations exist spanning different length scales, 
we may achieve a 'heavy tailed' statistics; that is, no 
convergence to a Gaussian and the distribution becomes 
Lorentzian. In reality this will always be a truncated 
Lorentzian due to the finite sample size, thus sidestep- 
ping any mathematical qualms regarding the use of the 
Lorentzian distribution. The Lorentzian is a scale inde- 
pendent distribution; consequently the number of steps, 
m, drops out of Eq. (|5c| . We note in passing that the vor- 
tex paths have a fractal character in this case. We have 
plotted the entropy for the uniform and angular functions 
in Appendix A. Some further details on the extrema of 
Eq. (|5bp are listed in Appendix B. 

This model naturally allows vortex interactions with 
multiple defect populations to be taken into account. 
For interactions with multiple Gaussian populations we 
have F 2 = T\ + F 2 , + . . ., while for multiple Lorentzians, 
r = Fi + F2 + . . .. These results are derived in Appendix 
C. Whether a particular peak is due to a unique popu- 
lation of defects or a compound population is a question 
of interpretation that needs to be addressed through the 
experimental evidence. 

As a method of inference, the principle of maximum 
entropy further allows us to make the following observa- 
tions. If we can describe J c (0) using Eqs. (0 then we 
have found the correct constraints for the problem and 
no better description of J c (0) is possible unless we can 
find weaker constraints and higher entropy functions to 
fit the data, which is highly unlikely. If J c (6) cannot be 
described by these equations then further constraints on 
the microstates exist which have not been taken into ac- 
count. The model thus provides a direct test of the high- 
field electron mass anisotropy scaling hypothesis, Eq. ([2]), 
since the model as presented takes no account of electron 
mass anisotropy. 



III. RESULTS AND DISCUSSION 

In this section we fit J c (9) data for a wide variety of 
defect-engineered samples, both LTS and HTS, to com- 
binations of Eqs. ([5]). In the vortex path model each 
equation corresponds to one set of constraints for one 
set of microstates. There is no expectation that a sin- 
gle equation enumerates all the possible microstates and 



so fitting involves choosing what, in the field of prob- 
ability, is termed a mixture distribution. For example, 
with both a6-plane and c-axis correlated defects present, 
both a so-called "a&-plane" peak and a "c-axis" peak in 
J c are possible. Thus, the combination of two orthogonal 
defect populations can lead to a mixture of two (or possi- 
bly more) distributions. The equations above have been 
written for a peak centered on the a&-plane; for a c-axis 
peak the 9 coordinate will be modified accordingly. Like- 
wise, although not explicitly considered here, it is per- 
fectly possible to account for inclined defect structures, 
for example due to off-axis irradiation, vicinal substrates, 
inclined substrate deposition or ion-beam assisted depo- 
sition, within the framework of the vortex path model 

The fitting has been done "by eye" with the assis- 
tance of a non-linear least squares fitting algorithm. In 
some data sets there exist asymmetries or drift which we 
have not attempted to correct. A formal statistical best 
fit procedure is possible but such a development is not 
warranted for the interpretation we are attempting here. 
More importantly we are attempting to interpret these 
distributions in terms of the underlying defect structure 
which contributes to each set of microstates or pinned 
vortex paths. This interpretation is of necessity tentative 
due to our limited knowledge of both the microstructure 
of the samples and the effectiveness of different elements 
of that microstructure in pinning. 



A. Nanostructured Nb films 

We begin by presenting data obtained on Nb thin films 
nanostructured to incorporate finely-tuned arrays of ver- 
tical pores. The sample preparation details and J C {B) 
data for these samples were published elsewhere [33| ■ The 
samples have a structure of nano-cnginecred pores etched 
vertically through the 36 nm thickness of the film with an 
average pore diameter of 75 nm and an inter-pore spacing 
of about 140 nm. The J c values for different samples are 
variable, but for the best samples approach the depairing 
current density Jd at self-field. The incorporation of the 
nano-pores increases J c by up to a factor of 50 over the 
virgin films, so clearly the pores contribute significant 
pinning. 

Measurements of J C (Q) on one of these samples are 
shown in Fig. Q] We note that this data is qualitatively 
similar to much of the literature data obtained on YBCO 
and observe that, following the approach of @ , such data 
would be interpreted as comprising three distinct angular 
regimes: (i) a sharp in-plane peak in the region of ±90°, 
attributed to staircase vortices predominantly following 
the intrinsic pinning a&-planes; (ii) a broader underlying 
in-plane peak attributed to mass anisotropy effects; and 
(iii) an enhanced J c near the out-of-plane direction at 
0° (in this case manifesting itself as a flat angular de- 
pendence over an extended angular range) , attributed to 
correlated out-of-planc pinning. 

Evidently, in the case of this data on Nb, intcrpreta- 
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FIG. 1. J C (P) at 5 K for a Nb thin film incorporating an array 
of vertical columnar pores: experiment (•), full fit (~~), fit 
components (— ,~ ) . The fitting parameters for these 
data are summarized in Table [T] 



TABLE I. Parameters of fit components in Fig. [T] 



Angular function 




0.1 T 


0.2 T 


Uniform (—) 


Jo 


1.140 




OOP Gaussian ( — ) 


Jo 

r 




0.351 
1.000 


IP Gaussian (— ) 


Jo 


0.015 


0.013 


r 


0.043 


0.030 


IP Gaussian (— ) 


Jo 


0.083 


0.426 


r 


0.249 


0.400 



tions (i) and (ii) are manifestly incorrect, and yet the 
features remain. Applying the vortex path model to this 
data, we also identify three components. However, our 
attribution of the origin of these components is somewhat 
different. In the vortex path model, each component is 
identified as resulting from the statistical combination 
of a dominant pinning population, giving correlated pin- 
ning aligned with the centre of the peak, and a pinning 
species providing a pinning force orthogonal to this di- 
rection, that broadens the peak. Thus, the two com- 
ponents (brown and green) centered on ±90° represent 
two distinct combinations of dominant in-plane pinning 
with a pinning population orthogonal to this (i.e. out- 
of-planc). Interestingly the two components have simi- 
lar peak height. We are not certain as to the origin of 
this effect but it may be due to the frequency of interac- 
tion of the vortex path with the orthogonal defects being 
similar. The two defect populations operate over differ- 
ent length scales (reflected in the widths of the peaks). 
Since this particular system is extremely simple, clean 
and well-defined, the only likely sources of correlated pin- 
ning are the surface interfaces (the film thickness being 



comparable to the coherence length) and the patterned 
vertical pores. We therefore identify the broader peak as 
arising from the interaction of the surface pinning with 
large-scale vertical defects comprising the pore structures 
themselves, while the narrow peak is the result of the in- 
teraction of the surface pinning with the (small-scale) 
surface roughness caused by the presence of the pores. 
It is then clear why these two interactions should occur 
with the same frequency. 

The third component in the low-field data appears in 
the form of a constant background. This is the uniform 
function. A better understanding of the origin of this 
component comes from considering the behavior under a 
larger applied field, where it is seen that this component 
narrows into a rather broad angular Gaussian. The angu- 
lar peak width of this component is similar to that of the 
broad in-plane peak, and so we attribute this third com- 
ponent to the dominant out-of-plane pinning facilitated 
by the pores, broadened by the interaction with the sur- 
face pinning. The scale factor is somewhat larger because 
instead of being limited by the length of the pore (i.e. the 
film thickness), it is now determined by the distance the 
vortex is able to travel in the in-plane direction while 
maintaining the macroscopic field orientation. Since this 
pinning mechanism is effective out to angles of about 
±60°, this distance is approximately tan 60° = -\/3 times 
the film thickness. At the lower field the weaker vortex- 
vortex interactions remove the statistical constraints on 
the majority of vortex paths formed between the pores 
and the surface pinning, thus resulting in the large uni- 
form contribution. 

It is to be emphasized from these results that the effect 
of the introduction of perfectly strong, perfectly corre- 
lated out-of-plane pinning defects is not necessarily the 
creation of a strong out-of-plane peak in the critical cur- 
rent. The absence of such a peak, as here, cannot there- 
fore be taken as evidence of the absence of correlated 
pinning. Here, the existence of such pinning centers has 
a determining effect on the critical current, but it does 
not manifest itself as a simple out-of-plane peak. Instead 
it contributes at all angles, through both the broaden- 
ing of the in-plane peak and the addition of a broad and 
almost flat component centered on the out-of-plane di- 
rection. This is of critical importance for the interpreta- 
tion of other data, for example in similar experiments to 
these recently conducted on YBCO [34[ which also lack 
an out-of-plane peak in J c (0). Likewise, 'shoulders' in 
the angular J c dependence at an intermediate angle be- 
tween the in-plane and out-of-plane directions are a nat- 
ural result of the interplay of in-plane and out-of-plane 
correlated pinning, and not a signature of some unortho- 
dox type of pinning at an arbitrary angle. As we shall 
see, such shoulders are an increasingly prominent feature 
of pinning-optimized samples, and one which the elec- 
tron mass anisotropy scaling approach is entirely unable 
to address. 
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FIG. 2. J c (0) at 77 K for a PLD YBCO film (data from [H): 
experiment (•), full fit (^), fit components (— ,—)■ The 
fitting parameters for these data are summarized in Table ITT1 



TABLE II. Parameters of fit components in Fig. [2] 



Angular function 




1 T 


3 T 


Uniform (— ) 


Jo 


0.600 


0.132 


ab Gaussian (—) 


Jo 


0.020 


0.004 


r 


0.075 


0.035 


ab Gaussian (— ) 


Jo 

r 


0.103 
0.473 


0.113 
0.396 



onal defects arc correlated with each other as was the 
case in the Nb films. We hypothesize that the explana- 
tion is the same: the population of longer c-axis defects 
(e.g. threading dislocations which penetrate the full film 
thickness) also creates small scale disorder (e.g. surface 
roughness) which vortices lying close to the a6-plane ori- 
entation see as a defect population available to provide 
the small amount of c-axis pinning required to broaden 
their peak. If we label the scale factors for these individ- 
ual populations as Ti and then the broad peak is most 
probably of scale factor T = (T^ + T\) 1 / 2 ^ i.e. we observe 
the contribution of pinning paths that interact with both 
species. Where the two distributions are strongly dis- 
parate, as here, the influence of the small scale disorder 
on the broader peak becomes negligible. Note that the c- 
axis defects present must of necessity also be contributing 
to the uniform component as this is the only component 
with a non-zero magnitude at 6 = 0. All other defects in 
the sample also contribute to the uniform component. 

At 3 T, the peak magnitudes of the narrow afe-peak and 
the uniform component are the same. The broad a&-peak 
is now of greater magnitude than these components. This 
is the same trend as was observed for Nb. This suggests 
that the only effective c-axis pinning remaining at this 
field is now provided by the through-thickness disloca- 
tions that match the frequency of the defects broadening 
the narrow a6-peak. These defects remain effective in 
broadening the afc-peak as at 1 T. However, the match- 
ing field of these defects is such that they are no longer 
effective pins in their own right at 3 T, and therefore the 
uniform component is diminished in magnitude. 



B. PLD YBCO films 

Turning now to HTS materials, we begin with the sim- 
plest data available for an epitaxial 300 nm thick YBCO 
film, without nanostructural modifications, prepared by 
pulsed laser deposition (PLD) on a single crystal SrTi0 3 
substrate [H|, and we see that all the features just de- 
scribed are again present (Fig. [2]). There is no need to 
invoke anything beyond the known microstructure of the 
sample to explain all of the observed features. In this 
thicker film, rather than surface pinning, effective in- 
plane pinning is provided by the intrinsic pinning due 
to the planar crystal structure of the sample, with non- 
superconducting planes separated by roughly the coher- 
ence length. Significant out-of-plane pinning is provided 
by the array of well-known structural defects resulting 
from the columnar growth mode of PLD samples that 
are the underlying cause of the high J c in thin films j36j . 

At 1 T, we observe three contributions to the pinning 
profile: two afr-centered angular Gaussians and a uni- 
form component. This indicates two distinct c-axis de- 
fect populations, each broadening the a6-peak by a dif- 
ferent amount. The peak heights of the two Gaussians 
again coincide, suggesting that the populations of orthog- 



C. MOD YBCO tapes 

The real power of the vortex path model, however, 
is revealed in its application to technologically-relevant 
HTS samples. In materials presently being developed for 
application, pinning effects dominate at all angles, neces- 
sitating consideration of the effects of combining multi- 
ple different defect populations. It is in these samples 
that the vortex path model comes into its own, explain- 
ing every aspect of the observed data in a coherent, self- 
consistent manner. 

Angular I c data for a production sample of American 
Superconductor Amperium™tape with a superconduct- 
ing YBCO layer fabricated by metal-organic deposition 
(MOD) are shown in Fig. [3] The major contributions 
to I c (&) can be fitted with three angular functions at all 
fields. A well-defined afe-plane peak increases in magni- 
tude as the applied field increases. This is paired with an 
extremely broad c-axis peak that broadens further under 
increasing field, resulting in the characteristic shoulder 
shape. Note that these shoulders are in approximately 
the same position as for the LTS film. The gradual broad- 
ening of the c-axis peak is consistent with the observation 
that the a6-plane pinning is increasing in strength, and 
will therefore contribute more to the broadening of the 
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FIG. 3. I c (0) at 77 K for MOD YBCO tape: experiment (•), 
full fit (—) , fit components (— ) . The fitting parame- 

ters for these data are summarized in Table 11111 Where there 
exists significant asymmetry in the data set, the first half of 
the data has been fitted. 



TABLE III. Parameters of fit components in Fig. [3] 



Angular function 




1 T 


2 T 


3 T 


Uniform (—) 


Jo 


117.8 


55.45 


28.15 


c Gaussian (— ) 


Jo 


0.331 


0.461 


1.755 


r 


0.056 


0.065 


0.133 


ab Gaussian (—) 


Jo 

r 


10.83 
0.664 


4.21 
0.962 


5.19 
1.69 



c-axis peak. Because the fit comprises just two angu- 
lar Gaussian functions centered on the orthogonal direc- 
tions, this indicates that there exists a single population 
of afr-plane defects and a single population of c-axis de- 
fects which are the dominant contributors to the angle- 
dependent pinning. The two peaks arise from the combi- 
nation of these same two pinning populations as outlined 
previously. 

The sample additionally has a significant uniform com- 
ponent at all fields. We have seen in the Nb data how this 
can arise from correlated defects. In these films it may 
also include a substantial contribution from point defects. 
An assembly of point defects may pin a single vortex or 
they may act in combination with correlated defects. To 
emphasize the importance of this interpretation, which 



differs from the standard view, we note that in the data 
measured at 1 T and 2 T there is a small angular region 
around ±75° where neither of the angular Gaussian func- 
tions has any magnitude. It is not credible to suppose 
that in precisely this region the correlated defects present 
in the sample provide no contribution to I Cl while con- 
tributing significantly at all other angles. As evidence to 
support this view, we have seen in the PLD sample that 
the afr-peak has a width extending over this region, while 
in this sample at 3 T we can again see directly the contri- 
bution of correlated defects in this region. However, we 
also know that in these pinning-engineered samples, the 
incorporation of nanoparticles — non-correlated defects 
— has increased the overall I c . Therefore both nanopar- 
ticles and correlated defects must be contributing to the 
uniform component. 

These samples differ from previously described MOD 
samples [26[ in having much stronger c-axis pinning. Pre- 
sumably this is due to twin planes or twin plane inter- 
sections, which are the only known high density c-axis 
defects to occur naturally in MOD samples. Again, the 
vortex path model enables us to relate the pinning pro- 
file observed in I c (0) to microstructural features of the 
sample that can be investigated experimentally. 

The fitting obtained for samples on technical sub- 
strates using angular Gaussians is seldom as good as 
for angular Lorcntzians, as seen throughout this paper 
and in (28|. This is because the samples are also disor- 
dered in their grain orientation, with an out-of-plane mo- 
saic spread around 5°. With the Lorentzian components, 
this actually helps to produce better fits as it generates 
the same result as having a collection of slightly differing 
length scales for the same defect population across dif- 
ferent grains. As explained in [28|, a Lorentzian is equiv- 
alent to a mixture of Gaussians and the disorder helps to 
create a spread of Gaussians. Conversely, when we have 
a Gaussian, the grain misorientation is working against 
a good fit. A Gaussian requires a single scale parame- 
ter but we don't have one because defects in differently 
oriented grains project onto our measurement axes differ- 
ently On this basis, it would be appropriate to smooth 
the fits over an angular range approximating the mosaic 
spread in the sample. Further counting against a good 
Gaussian fitting, we will never truly have a single defect 
population giving all the correlated J c (0) — there are 
bound to be minor defect populations also contributing. 
This is most likely the case in this sample where in the 2 
T and 3 T data a second c-axis peak can just about be 
resolved. Nonetheless, even without accounting for these 
finer details, we obtain a high level of accuracy with very 
few components. 



D. YBCO with Ba 2 YNb0 6 additions 

We now turn to defect-engineered samples in the true 
meaning of the term; that is to say samples with incor- 
porated defects designed to influence the pinning pro- 
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FIG. 4. J c (0) at 77 K for a PLD YBCO film with Ba 2 YNb0 6 
additions (data from [35|]): experiment (•), full fit (—), fit 
components (— ,~ ,~ ) . The fitting parameters for these 
data are summarized in Table [TV] Where there exists signifi- 
cant asymmetry in the data set, the first half of the data has 
been fitted. 



TABLE IV. Parameters of fit components in Fig. 2] 



Angular function 




1 T 


3 T 


Uniform (—) 


Jo 


1.040 


0.273 


c Lorentzian / Gaussian {——) 


Jo 

r 


L 0.223 
L 1.012 


G 0.075 
G 0.524 


ab Lorentzian / Gaussian (—) 


Jo 

r 


L 0.050 
L 0.302 


G 0.138 
G 0.050 


ab Gaussian (—) 


Jo 

r 




0.125 
0.570 



file in particular ways. The first of these we will con- 
sider is the introduction of self-assembled c-axis aligned 
nanocolumns of second-phase material intended to pro- 
vide strong, effective flux pinning in the c direction where 
J c is usually lowest. In the present case, this has been 
achieved through the addition of the double-perovskitc 
Ba 2 YNb0 6 to PLD YBCO films Q. This material is 
incorporated into the YBCO matrix as ^100 nm long, 
10-15 nm diameter rods aligned with the c-axis of the 
YBCO and has enabled a two to three times enhance- 
ment in J c (B\\c) compared with pure YBCO films. The 
data shown in Fig. 0] arc for a 300 nm thick film compris- 
ing 5 mol% Ba2YNbOg, resulting in nanocolumns having 
an average lateral spacing of around 40 nm. 

At low field (1 T), we once again require three com- 
ponents to fit the data: a uniform component and two 
Lorentzian peaks. The presence of Lorentzian rather 
than Gaussian peaks may be caused by additional lattice 
disorder introduced into the YBCO by the nanocolumns 
due to their large lattice mismatch and consequent strain- 



ing of the matrix. Significant distortion of the YBCO 
lattice has been observed in samples of this composition. 
Additionally, splaying of the nanorods about the c-axis 
means the direction of their pinning is not well-defined, 
but instead encompasses a distribution of its own. Con- 
sequently, instead of a single defect population yielding 
a single scale parameter there are multiple Gaussians 
present which have merged to produce a Lorentzian line- 
shape. In contrast to the pure YBCO films, here the 
two peaks are not both afr-plane peaks, but rather one 
of them is a strong c-axis peak as a result of the highly 
effective nanocolumnar pins. In fact, at this field, these 
pins are so effective that the magnitude of the c-axis peak 
is greater than that of the a&-plane peak. This in turn 
explains the absence of the second a6-plane peak in this 
data: the strength of the c-axis pinning is such that it is 
completely dominant in its interaction with the available 
a6-pinning. This results in the strong uniform compo- 
nent ranging across the angles usually occupied by the 
broad et6-plane peak. 

At higher field (3 T), the nanocolumnar pinning begins 
to lose its strength, a common phenomenon in samples 
utilizing present flux pinning strategies. This is because 
the overall density and therefore matching field is rela- 
tively low for this defect population. Clearly the pin- 
ning centers are still present, but they lose in dominance 
to the afc-planc pinning, and the broad peak associated 
with their interaction with this pinning source is recov- 
ered. The response at this field is extremely similar to 
the pure YBCO PLD films, with the addition of a small 
c-axis peak due to the residual effect of the additional 
engineered c-axis pinning. The higher field thus again 
helps to distinguish the different pinning contributions. 
The change in peak shape from Lorentzian to Gaussian 
likely results from vortex- vortex interactions at high field 
acting to constrain the vortices, diminishing both the ef- 
fectiveness of the pinning and reducing the effect of splay. 
We also see the previously noted phenomenon of the two 
a6-plane peaks coinciding in magnitude, and additionally 
here the uniform component also coincides. The uniform 
component is about twice the magnitude seen in the pure 
YBCO, reflecting the continuing contribution of the en- 
gineered c-axis pinning to the uniform component, not 
just to the c-axis peak. 



E. YBCO with Gd 3 Ta0 7 additions 

In the attempt to progress beyond BaZrOs as an ef- 
fective second-phase pinning inclusion for YBCO, and to 
increase the effectiveness of this type of artificial pinning 
center, the pyrochlore rare earth tantalate i?3Ta07 ma- 
terial was successfully introduced, resulting in significant 
J c enhancement compared to BaZrOa-containing sam- 
ples across the angular range [38$ . Due to its close lattice 
match to YBCO, in contrast to BaZr03 and Ba2YNbOe, 
i?3Ta07 addition readily produces extremely fine and 
highly linear through-thickness nanorods without splay- 
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ing, and with no poisoning effect reducing the critical 
temperature. By a simple variation in deposition rate, 
the nanostructure of the i?3Ta07 additions can be con- 
trolled, either freezing the additions into a fine dispersion 
of nanoparticlcs of ~5 nm diameter (at high deposition 
rate) , or allowing time for diffusion-mediated growth into 
extended nanorods of 5-10 nm diameter and 10-15 nm 
spacing (at low deposition rate). This provides a unique 
'lever' offering unprecedented control over the nature of 
the flux pinning in samples incorporating this additive 

The data shown in Fig.[S]for 500 nm thick YBCO films 
prepared by PLD with 1.5 mol% GdaTaOy additions ex- 
emplify the pinning control that can be achieved with this 
additive, where the anisotropy is seen to vary from broad 
afr-peak in the case of a high-rate deposition through to 
strong c-axis peak for low-rate deposition. Supported 
by TEM imaging, the simple additive model of pinning 
used in [39| to describe this data provides an adequate 
explanation of the primary features, that is to say the ex- 
istence of a broad, flat profile with residual ab-peaks in 
the case of the isotropic nanoparticlc pinning produced 
at high growth rate giving way to a strong c-axis peak in 
the presence of the nanocolumns formed at low growth 
rate. However, application of the vortex path model to 
these data allows us to decompose the behavior further. 

In the case of the sample deposited at high rate, we see 
the simplest behavior observed to date — just two com- 
ponents: one uniform and one broad et6-peak. Similar 
to the pure PLD YBCO film, we attribute the broad ab- 
peak to dominant intrinsic a6-planc pinning broadened 
by the type of through-thickness defects characteristic 
of PLD. The effect of the nanoparticlc inclusions is to 
strongly enhance the uniform component with the result 
that the sharp afr-peak previously observed in films of this 
type is obscured. This is an example of similar-strength 
pinning centers adding non-linearly to generate the total 
pinning force, and evidences the failure of a simple ad- 
ditive pinning model. As the applied field is increased, 
the effectiveness of the nanoparticulatc pinning centers 
diminishes, leading to a reduction in the uniform compo- 
nent, and the sharp afr-plane peak is again distinguished. 

Conversely, in the case of the sample deposited at low 
rate, the most prominent feature is the c-axis peak result- 
ing from the nanorods. The vortex path model represents 
this as three angular Gaussians of increasing width, all 
centered on the c-axis. Here, due to the extreme linearity 
and crystalline perfection of the nanorods, no tendency 
towards a Lorentzian lineshape is observed, and distinct 
broadening defect populations can be resolved. Three 
components of differing width imply three distinct pop- 
ulations of broadening defects available to interact with 
the nanorods. Based on the limited microstructural data 
available on these samples, we speculate that, in common 
with the MOD samples, strong a6-plane pinning is one 
of the broadening populations, while another broaden- 
ing defect may be specific a&-plane defects such as stack- 
ing faults. These tentative assignments can be verified 



or refuted through experiments designed to control the 
density of specific defects within the sample and correlate 
this with J c (0). As the field is increased, these three com- 
ponents become better resolved, broadening and devel- 
oping the characteristic shoulders also seen in the MOD 
samples as the intrinsic afe-plane pinning becomes rela- 
tively more effective than the artificial c-axis pinning. 

The advantage of the vortex path model in analyzing 
this data is that we can see how the basic components of 
correlated and uncorrelated defects are combining to pro- 
duce the overall J c (9) behavior. That particular peaks 
are narrow or broad does not mean that whole pinning 
populations are suddenly appearing or disappearing but 
rather that the relative strengths and the statistics of 
how these populations combine is being altered. 



F. YBCO with Ba 2 YNbO e + Gd 3 Ta0 7 additions 

Recently, the importance to pinning optimization of 
achieving segmented nanocolumnar pins rather than con- 
tinuous nanocolumns was suggested floj to mitigate the 
effects of thermally activated depinning. One method 
by which this has been achieved is through the simul- 
taneous incorporation of the two types of pinning ad- 
ditive already considered, Ba2YNbOg and GdaTaOy. 
What forms in this case is a complex composite of 
YBCO, Ba2Y(Nb,Ta)06 segmented nanorods, in-plane 
Y2O3 platelets and nanoparticles of the YBa2Cu40g su- 
perconducting phase, where the Gd partially substitutes 
for Y in all four phases and the Ta and Nb content goes 
entirely into forming nanorods [35j. The J c {9) depen- 
dence as a function of the applied field for a 300 nm thick 
sample of this composition prepared by PLD is shown in 
Fig.© 

From the low field data, the most striking feature is 
the large sharp c-axis peak, far stronger and better de- 
fined than we have seen in any other data. Analysis in 
terms of the vortex path model begins to give us some 
idea of why this is — the coincidence of no less than 
four components centered on the c direction combine to 
give this peak. Its sharpness is caused by the combi- 
nation of two narrow components of similar magnitude, 
while its baseline is boosted by the contribution of two 
additional broad components, as well as a large uniform 
component. The effect of this large uniform component 
is to make the a&-peak look small, but once again this 
is not due to any particular lack of a&-pinning in these 
samples, but rather the strong relative dominance of the 
c-axis pinning. In the same way that strong a6-pinning 
can obscure a c-axis peak despite the existence of signif- 
icant c-axis pinning, uncommonly strong c-axis pinning 
as has been achieved here can also obscure the afr-peak. 
The uniform component is similar in magnitude to the 
nanorod-containing samples of the previous section. 

By analogy with the earlier niobate addition data, we 
identify the broader of the two sharp c-axis peaks as re- 
sulting from the segmented nanorods, in interaction with 
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FIG. 5. J c (6) at 77 K for PLD YBCO films with GdsTaOy additions, deposited at different rates to induce nanoparticle or 
nanorod formation (data from [13] ): experiment (•), full fit (—), fit components {^—,^—,^—,^—,^—,—)- The fitting parameters 
for these data are summarized in Table IVl 



TABLE V. Parameters of fit components in Fig. [5] 



Angular function 




Hit 
1 T 


;h rate 
3 T 


Low rate 
IT 3 T 


Uniform (—) 


Jo 


2.293 


0.534 


1.084 


0.669 


c Gaussian (— — ) 


Jo 

r 






0.117 
0.263 


0.043 
0.244 


c Gaussian (— — ) 


Jo 






0.207 


0.055 


r 






0.630 


1.100 


c Gaussian ( ) 


Jo 

r 






0.099 
1.444 


0.027 
2.816 


ab Gaussian (—) 


Jo 

r 




0.050 
0.154 


0.011 
0.045 


0.003 
0.030 


ab Gaussian (—) 


Jo 


0.253 


0.128 






r 


0.384 


0.465 







the intrinsic a&-plane pinning, observing that it varies 
with field in the same way, diminishing and broadening 
out to a relatively flat response by 3 T. This leaves the 
sharper c-axis peak to be identified, and we note that it 
tracks the broader peak in magnitude. It also remains rel- 
atively constant throughout the entire field range, chang- 
ing little in width or relative magnitude. This suggests 
two things: firstly that it has a high matching field, i.e. 
it results from a high-density field-independent popu- 
lation of defects, and secondly that its width probably 
arises from a property such as the angular spread of the 
nanorods. Hence this is an unusual peak in that it may 
arise from a single defect species only, interacting with 
itself through its own disorder. 

The broader c-axis components then likely arise due 
to interactions between the segmented nanorods and the 
in-plane Y2O3 platelets and YBa 2 Cu40g nanoparticles. 
These unique features are a direct result of the unique 



microstructure of these samples, with the Y2O3 platelets 
being of similar dimension to the nanorod segments. It 
is therefore unsurprising that they greatly extend the an- 
gular range over which the nanorods can effectively pin, 
facilitating pinning of vortices between adjacent nanorod 
segments at high field angles. 

These four c-axis components can be traced through to 
3 T, by which point the c-axis pinning is significantly re- 
duced in strength, as is commonly observed for all artifi- 
cial pinning centers known to date. The two broad c-axis 
components continue broadening strongly with increas- 
ing applied field. To higher fields, the dominance of the 
c-axis pinning is overcome, and familiar features relating 
to the a&-pinning begin to emerge. The resultant combi- 
nation of these similar strength features yields extremely 
complex angular variations in J c that arc nonetheless re- 
solvable into a few primitive components that describe 
the subtlest features of the data extremely well. 

At 4 T, an interesting crossover occurs. It is at this 
field that the strongest nanorod component is broadened 
to such an extent that it has developed significant shoul- 
ders that result in a quite notable response in the J c {&)- 
It is appropriate that such significant microstructural fea- 
tures should continue to give a marked response in the 
J c (9) even though their contribution to the overall J c is 
now diminished. At the same time, the broadest c-axis 
peak vanishes, to be replaced with the usual afe-pinning 
peak that we attribute to intrinsic pinning being broad- 
ened by through-thickness defects. The disappearance is 
assumed to be a consequence of the applied field becom- 
ing greater than the matching field of the relevant c-axis 
defect species. 

We note a couple of interesting points about the fit- 
ting. Again we observe coincidences in the peak heights 
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FIG. 6. Jc(6) at 77 K for a PLD YBCO film with Ba 2 YNb0 6 + Gd 3 Ta0 7 additions (data from [H; 2 T data unpublished 
data from the same sample): experiment (•), full fit (— ), fit components (^— ) . The fitting parameters for 
these data are summarized in Table [VTl The small afo-peak has been ignored in the fitting. 



TABLE VI. Parameters of fit components in Fig. [6] 



Angular function 




1 T 


2 T 


3 T 


4 T 


5 T 


6 T 


Uniform 




Jo 


1.166 


0.829 


0.801 


0.531 


0.311 


0.148 


c Gaussian 




•A. 


0.125 


0.105 


0.070 


0.033 


0.035 


0.025 




r 


0.091 


0.113 


0.176 


0.153 


0.151 


0.152 


c Gaussian 




Jo 


0.348 


0.384 


0.356 


0.256 


0.185 


0.118 




r 


0.257 


0.371 


0.722 


1.21 


1.43 


1.55 


c Gaussian 




Jo 


0.404 


0.368 


0.163 


0.045 


0.023 


0.015 




r 


0.672 


0.831 


1.49 


3.56 


4.00 


3.64 


c Gaussian 




Jo 


0.218 


0.140 


0.055 










r 


1.63 


2.38 


3.50 








ab Gaussian 




Jo 








0.070 


0.118 


0.145 




r 








0.491 


0.422 


0.400 



of multiple components. These are seen to occur not only 
at the center of the distribution but also at other points 
where the distributions reach a local extrcmum. At an 
cxtremum dg/d6 = and the distribution is locally un- 
constrained. This coincidence may therefore have an in- 
terpretation in terms of the maximization of entropy for 
a mixture distribution which we do not presently under- 
stand. The other point we make is to remind the reader 
that each component has only a single scale parameter. 
This means that both the shape and the position of the 
shoulders are highly constrained by the form of the angu- 
lar Gaussian distribution we have derived. It is difficult 



to conceive of other primitive functions which could so 
well reproduce these features of this data. 

IV. CONCLUSION 

We have addressed the question of whether the applica- 
tion of the electron mass anisotropy scaling approach to 
field angular J c {6) data of pinning-optimized samples is 
theoretically justified, and have concluded that it is not. 
In experimental support of this determination, we have 
presented data measured on thin films of the isotropic 
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superconductor Nb that exhibit features commonly at- 
tributed to electron mass anisotropy, reinforcing the fact 
that J c is a wholly extrinsic material parameter and de- 
manding an alternative explanation that can be related 
to the microstructure or other extrinsic properties of the 
sample. We provide this explanation in terms of the vor- 
tex path model, and have applied the same model to a 
wide range of contemporary pinning-engineered samples, 
from which a number of general results emerge. Firstly, 
that the introduction of perfectly strong, perfectly corre- 
lated out-of-planc pinning defects docs not necessarily re- 
sult in a strong out-of-plane peak in the J c {9), and there- 
fore that the absence of such a peak cannot be taken as 
evidence of the absence of correlated pinning. Secondly, 
that pinning from different sources is not a simple addi- 
tive process, but rather depends on what pinning is al- 
ready present. In this context, common manifestations of 
strong out-of-plane correlated pinning are a broadening 
of the in-plane peak or an almost flat angular variation in 
J c (0) about the out-of-plane direction, while commonly- 
observed shoulders in the angular J c dependence are a 
result of the interplay between in-plane and out-of-plane 
correlated pinning of similar strength, and not a signa- 
ture of exotic pinning at an arbitrary angle. In particular 
we emphasize that a broad in-plane peak is evidence of 
strong out-of-planc correlated pinning, and not a reflec- 
tion of electron mass anisotropy. 

Our assignment of defect populations to particular 
peaks is based on the available knowledge of the sample 
microstructure. Correct identification of the individual 
pinning populations requires a correlation between J c (9) 
and microstructural analysis over a family of samples in 
which the density of these defect populations is system- 
atically varied and the effect on the different components 
of J c quantified. These hypotheses can then be tested by 
measuring samples in those temperature and field ranges 
which best distinguish between the different defect pop- 
ulations. It appears from our analysis of available data 
that measuring the field dependence of the angular profile 
is very useful for determining the origin of pinning effects. 
The evolution of the magnitude of a component allows us 
to determine the matching field for a contributing defect 
population. The evolution of the shape also correlates 
with matching field effects, e.g. out-of-plane components 
broaden towards the in-plane direction as the field in- 
creases due to the high matching field of the intrinsic 
pinning. 

In the context of the principle of maximum entropy, 
the success of the vortex path model in describing the 
data shows that no further constraints beyond those we 
have imposed on the vortex path statistics are necessary 
to explain J c (9). Thus we can conclude with confidence 
that the electronic mass anisotropy plays no role in deter- 
mining the form of J c {9) for any of the various samples 
studied. If we find samples which do not immediately 
fit the vortex path model equations then different con- 
straints are operating in those systems. An obvious ex- 
ample is the need to extend the model to cases where 



J c {9) has asymmetries due to oblique correlated defect 
populations. The success of the vortex path model in 
this instance also suggests that other angular data in- 
fluenced by the pinning landscape could be better de- 
scribed by Eqs. (O, or at least other expressions derived 
from maximum entropy considerations. The Ginzburg- 
Landau mass anisotropy expression is often applied sim- 
ply because no alternative angular expressions have been 
proposed. For example, we can see no reason why the 
angular dependence of the irreversibility field or the flux 
flow resistivity should follow an expression involving elec- 
tronic mass anisotropy. 
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Appendix A: Entropy of the angular distributions 

The entropy of the normalized uniform, angular Gaus- 
sian and angular Lorentzian distributions derived as 
Eqs. ([5]) is calculated from Eq. (|4|). This is done nu- 
merically as no analytical solution exists. For discrete 
distributions H > but for continuous distributions this 
restriction is relaxed with the understanding that all en- 
tropies are relative quantities. The results arc plotted 
on Fig. [7] Most striking is the special case whereby the 
entropy of the angular Lorentzian with scale factor T = 1 
equals that of the uniform distribution — a set of vor- 
tices constrained to pin along correlated defects can, in 
principle, result in a critical current distribution having 
the same entropy as that resulting from a set of fully 
unconstrained vortices. 

Appendix B: The extrema of Eq. (|5b|) 

Equation (|5b|) has extrema at 6 = and 9 = ±7r/2, 
and for T > 1/V2, at 9 = sin -1 (1/V2T). The magnitude 
at 9 = is J c = 0, the global minimum. At 6 = ±ir/2, 
the magnitude is J c = Jq/T, which is a maximum for 
r < l/\/2 and a local minimum for T > l/\/2. At 9 = 
sin-^l/v^r), J c = 2J rexp((l/2r 2 ) - 1), a maximum. 

Appendix C: Summation of pinning populations 

Suppose we have multiple pinning populations broad- 
ening a peak which individually in the vortex path model 
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Scale factor r 

FIG. 7. Entropy of the uniform, angular Lorentzian and angu- 
lar Gaussian distributions in dependence on their scale factor. 
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